APPROVED  FOR  PUBLIC 


Prepared  by 
Helmut  Y.  A.  Hsiao 


TEXAS  INSTRUMENTS  INCORPORATED 
Equipment  Group 
Post  Office  Box  6015 
Dallas,  Texas  75222 


Prepared  for  ^ 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
Alexandria,  Virginia  22314 


Sponsored  by 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
Nuclear  Monitoring  Research  Office 
ARPA  Program  Code  No.  7F10 
ARPA  Order  No.  2551 


Acknowledgment:  This  research  was  supported  by  the  Advanced 
Research  Projects  Agency,  Nuclear  Monitoring  Research  Office, 
under  Project  VELA- UNIFORM,  and  accomplished  under  the  tech 
nical  direction  of  the  Air  Force  Technical  Applications  Center 
under  Contract  Number  F08606-77-C-0004. 


0 4 02  0] 


Equipment  Group 

0 


V- 


APPLICATION  OF  A COMBINED  SOURCE  MODEL  FOR  SEISMIC  DISCRIMINATION 


TECHNICAL  REPORT  NO.  21 

VELA  NETWORK  EVALUATION  AND  AUTOMATIC  PROCESSING  RESEARCH 

Prepared  by 
Helmut  Y.  A.  Hsiao 


TEXAS  INSTRUMENTS  INCORPORATED 
Equipment  Group 
Post  Office  Box  6015 
Dallas,  Texas  75222 

Prepared  for 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
Alexandria,  Virginia  22314 

Sponsored  by 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
Nuclear  Monitoring  Research  Office 
ARPA  Program  Code  No.  7F10 
ARPA  Order  No.  2551 


25  October  1978 


Acknowledgment:  This  research  was  supported  by  the  Advanced 
Research  Projects  Agency,  Nuclear  Monitoring  Research  Office, 
under  Project  VELA- UNIFORM,  and  accomplished  under  the  tech- 
nical direction  of  the  Air  Force  Technical  Applications  Center 
under  Contract  Number  F08606-77-C-0004. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAPE  (Whmr  Data  EmaraQ 

I REPORT  DOCUMENTATION  PAGE 


a.  OOVT  ACCESSION  NOJ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
s.  recipient's  catalog  number 


S.  TYNE  OF  REPORT  A PERIOD  COVERED 


^APPLICATION  OF  A COMBINED  JIOURCE 
"MODEL  FOR  SEISMIC  DISCRIMINATION  „ 


<?v  Technical 


r /.,  i fiA, 


, -ALEX  (01 1-  TR  - 78 -J09 


17.  author^ 


(Helmut  Y.  A.  /ti siao  ! 


F08606-77-C-0004 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  • WORK  UNIT  NUMBERS 


V PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  10  PROGRAM  ELEMENT.  PRC 

Texas  Instruments  Incorporated  * wom<  UH,T  NUM 

Equipment  Group  / VELA  T/8705/B/PMP 

Dallas,  Texas 75222 

II.  CONTROLLING  OFFICE  NAME  AND  ADORESS  S **.  REPORT  PATE 

Advanced  Research  Projects  Agency  f 1 L f?§~Octi»5wSW7l 

Nuclear  Monitoring  Research  Office  \ A UUIHIR  5P  PISli 

Arlington,  Virginia 22209 101 i 

14.  MONITORING  AGENCY  NAME  S AOORSSSfU  dUlaNnl  from  Controlling  Olllem)  IS.  SECURITY  CLASS,  (ol  Ihfl 


(3* 


lAz. 


Air  Force  Technical  Applications  Center 
VELA  Seismological  Center 
Alexandria,  Virginia  22314 


UNCLASSIFIED 
Sa.  declassification/downgraoino 

SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  ami*  RaparO 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


17.  DISTRIBUTION  STATEMENT  (•!  Ma  aAalrael  anlara*  In  Slant  10.  II  RMWawl  fcaai  Wafa, 

(Tl4  rv  - -77- a ■ - ^ / 


I IS.  supplementary  notes 


ARPA  Order  No.  2551 


■ Ac  . -tsl  ~-i- 


It.  KEY  WORDS  (Carllnut  an  ravaeaa  alpa  II  naeaaaaiy  anS  MaiUiy  t y Slant  numht) 

Long-period  teleseismic  surface  wave  Moment  versus  plot 
Seismic  source  parameter  estimation  Amplitude  spectral  fitting 
Seismic  source  mechanism  Seismic  discrimination 

.^Combined  source 

10.  A»iTfr»qT  (Contlnu*  an  f >lrf»  If  wcsusry  m4  iSZmtitr  Mock  iwwfcfj 

^The  purpose  of  this  study  has  been  to  use  the  observed  long-period 
teleseismic  surface  wave  data  of  earthquakes  and  presumed  underground 
explosions  to  determine  the  source  mechanisms  and  to  use  these  source 
mechanism  estimates  to  discriminate  between  the  two  classes  of  events. 

In  order  to  demonstrate  this  purpose,  fifteen  events  have  been  selected  and 
analyzed  for  the  Russian  Eastern  Kazakh  (EKZ)  region,  central  and 
northern  Eurasia,  and  the  United  States  Nevada  Test  Site  (NTS).  — > 

DO  •JAN>7l  1473  EDITION  OF  I NOV  •»  IS  OBSOLETE  UNCLASSIFIED 


07Cr 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  Data  EmaratfJ 


UNCLASSIFIED 

MCUWITV  CLAMiriOTtOW  or  thi»  namto  balm  ImiwO 


C. 61 


20.  continued 

> The  observation  stations  available  for  this  study  included  the  Seismic 
Research  Observatories  (SRO),  Abbreviated  Seismic  Research  Observatories 
(ASRO),  the  High  Cain  Long- Period  (HCLP)  stations  (the  same  as  the  pre- 
vious VLPE),  and  existing  VELANET  arrays. 


For  a fair  comparison  of  explosion  and  earthquake  events,  both  have 
been  treated  in  the  same  way  by  using  the  combined  source  model  which 
combines  a point  explosive  source  and  a point  double-couple  source. 

c In  order  to  cut  down  computer  expenses  and  widen  the  searching  range 
of  some  parameters,  a binary  exhaustive  search  method  was  applied  to 
modify  the  amplitude  spectral  fitting  process  without  any  constraint.  The  re- 
sults obtained  from  this  modified  process  gave  very  good  spectral  fits  for 
most  events  and  lowered  the  minimum  residual  fitting  error  of  each  event. 

The  estimated  seismic  moment  versus  m^  plot  has  been  found  to  be 
very  promising  for  the  purpose  of  discrimination  between  earthquakes  and 
explosions.  Taking  the  surface  wave  seismic  moment  to  be  a universal 
scaling  parameter  for  sizing  events,  the  average  m^  of  explosions  is  one 
magnitude  greater  than  the  average  m^  of  earthquakes.  As  expected,  this 
indicates  more  efficient  radiation  of  high  frequency  P-waves  for  explosions 
of  equal  moment.  Since  this  criterion  appears  to  apply  equally  to  both 
Nevada  Test  Site  and  Eastern  Kazakh  presumed  explosions,  it  suggests  that 
surface  wave  moment  estimates  may  be  a more  universal  parameter  than  Mg 
for  scaling  the  size  of  earthquakes  and  explosions.  This  inference  from 
our  results  needs  further  investigation. 


UNCLASSIFIED 


tCCuMirr  CkMMPtCSTMM  Of  THIS  SftSSflntM  Dm,  Baton*) 


ABSTRACT 


The  purpose  of  this  study  has  been  to  use  the  observed  long- 
period  teleseismic  surface  wave  data  of  earthquakes  and  presumed  underground 
explosions  to  determine  the  source  mechanisms  and  to  use  these  source  mech- 
anism estimates  to  discriminate  between  the  two  classes  of  events.  In  order  to 
demonstrate  this  purpose,  fifteen  events  have  been  selected  and  analyzed  for  the 
Russian  Eastern  Kazakh  (EKZ)  region,  central  and  northern  Eurasia,  and  the 
United  States  Nevada  Test  Site  (NTS). 

The  observation  stations  available  for  this  study  included  the  Seis- 
mic Research  Observatories  (SRO),  Abbreviated  Seismic  Research  Observa- 
tories (ASRO),  the  High  Gain  Long- Period  (HGLP)  stations  (the  same  as  the 
previous  VLPE),  and  existing  VELANET  arrays. 

For  a fair  comparison  of  explosion  and  earthquake  events,  both 
have  been  treated  in  the  same  way  by  using  the  combined  source  model  which 
combines  a point  explosive  source  and  a point  double-couple  source. 

In  order  to  cut  down  computer  expenses  and  widen  the  searching 
range  of  some  parameters,  a binary  exhaustive  search  method  was  applied  to 
modify  the  amplitude  spectral  fitting  process  without  any  constraint.  The  re- 
sults obtained  from  this  modified  process  gave  very  good  spectral  fits  for  most 
events  and  lowered  the  minimum  residual  fitting  error  of  each  event. 

The  estimated  seismic  moment  versus  m^  plot  has  been  found 
to  be  very  promising  for  the  purpose  of  discrimination  between  earthquakes 
and  explosions.  Taking  the  surface  wave  seismic  moment  to  be  a universal 
scaling  parameter  for  sizing  events,  the  average  m^  of  explosions  is  one 
magnitude  greater  than  the  average  of  earthquakes.  As  expected,  this 


indicates  more  efficient  radiation  of  high  frequency  P-waves  for  explosions 
of  equal  moment.  Since  this  criterion  appears  to  apply  equally  to  both  Nevada 
Test  Site  and  Eastern  Kazakh  presumed  expkmons,  it  suggests  that  surface 
wave  moment  estimates  may  be  a more  universal  parameter  than  M for  scal- 
ing the  size  of  earthquakes  and  explosions.  This  inference  from  our  results 
needs  further  investigation. 
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SECTION  I 
INTRODUCTION 


The  primary  intention  of  this  report  is  to  obtain  a better  under- 
standing of  earthquake  and  explosion  source  mechanisms  as  determined  from 
long-period  surface  wave  spectra  and  to  apply  these  results  to  discriminating 
between  earthquakes  and  explosions.  Tsai's  method,  an  exhaustive  comparison 
scheme  between  the  observed  and  theoretical  spectra,  was  used  and  modified 
by  applying  the  binary  method  to  search  through  the  following  three  source 
parameters;  depth,  strike  angle,  and  F value  (F  is  the  ratio  of  the  point  radial 
and  double-couple  source  terms).  The  purpose  of  this  modification  is  to  econo- 
mize computer  expenses  and  to  acquire  more  precise  source  parameter  esti- 
mates over  a wider  range  of  the  F value.  In  order  to  improve  the  source  rep- 
resentation, a combined  source  model  which  consists  of  a point  double-couple 
source  and  a point  explosive  source  was  utilized  to  analyze  all  the  events  se- 
lected. The  following  report  represents  the  results  of  study  on  long-period 
surface  waves  of  fifteen  EKZ,  Eurasian,  and  NTS  events. 

Section  II  offers  a detailed  description  of  the  data  selection,  data 
processing,  data  availability  and  data  quality  of  this  study.  In  Section  HI,  a 
brief  summary  of  the  theoretical  source  description,  the  procedures  to  obtain 
observed  surface  wave  spectra,  and  the  amplitude  spectra  fitting  process  is 
presented.  The  results  of  source  mechanism  and  source  discrimination  are 
discussed  in  Section  IV.  Finally,  in  Section  V,  the  directions  this  study  expects 
to  take  in  the  future  are  discussed. 


SECTION  II 
DATA  BASE 


A.  DATA  SELECTION 

In  the  previous  report  (Sun,  1977),  twenty-seven  underground 
nuclear  explosions  from  three  different  regions,  the  United  States  Nevada 
Test  Site  (NTS),  the  Russian  eastern  Kazahk  (EKZ),  and  North  Caspian  Sea 
area,  were  selected  by  means  of  chronologically  exhaustive  search  (starting 
from  1976  and  progressing  backward  into  1975,  1974  and  so  on)  through  the 
Norwegian  Seismic  Array  (NORSAR)  bulletin  and  Preliminary  Determination 
of  Epicenters  (PDE)  bulletin.  The  availability  and  quality  of  NTS  events  data 
are  much  better  than  that  of  the  EKZ  events  data.  Therefore,  in  this  year's 
study,  only  EKZ  events  as  a group  of  presumed  explosions  were  searched 
(covering  1976  to  1977)  using  the  National  Earthquake  Information  Service 
(NEIS)  list  and  PDE  bulletin.  In  the  same  way,  an  attempt  was  made  to  search 
and  select  as  many  recent  Eurasian  earthquakes  as  possible  to  satisfy  the  pur- 
pose of  this  study. 

In  the  preliminary  list,  there  were  eighteen  EKZ  events  (eight 
in  1976  and  ten  in  1977)  and  twenty-one  Eurasian  events  (nine  in  1976  and  twelve 
in  1977).  However,  considerable  numbers  of  those  events  were  lost  due  to 
operational  and  recording  problems,  the  details  of  which  are  discussed  in 
Subsection  C.  Therefore,  two  processed  earthquakes  and  five  presumed  und- 
erground nuclear  explosions  were  selected  from  previous  reports  (Lambert  et 
al. , 1974  and  Sun,  1977)  to  supplement  the  data  base.  Based  on  these  process- 
ed events  and  the  preliminary  list,  the  availability  of  good  quality  long-period 
teleseismic  surface  wave  data  was  examined  and  the  final  list  of  the  selected 
events  was  obtained.  This  final  list  of  events  is  depicted  in  Table  U-l. 
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LIST  OF  THE  SELECTED  EVENTS  AND  THEIR  DATA  AVAILABILITY 
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Both  Rayleigh  and  Love  Wave  R:  Only  Rayleigh  Wave  L:  Only  Love  Wave 
Event  Selected  from  Previous  Report  (Use  the  same  event  names  as  before) 


LIST  OF  THE  SELECTED  EVENTS  AND  THEIR  DATA  AVAILABILITY 
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B. 


DATA  PROCESSING 


After  selecting  the  events,  the  next  effort  was  to  obtain  the  de- 
sired surface  wave  amplitude  spectra  from  ten  to  fifty  seconds  for  the  estima- 
tion of  the  seismic  source  parameters.  In  order  to  obtain  these  amplitude 
spectra,  the  field  tapes  which  contain  the  surface  wave  recordings  for  the  se- 
lected events  at  all  available  observation  stations  were  processed  as  follows: 

• The  desired  surface  wave  recordings  on  the  field  tape  were 
edited  to  yield  the  event  signal  which  usually  consists  of  3-N 
signal  traces,  where  N is  the  number  of  recording  sites  in  a 
station,  i.  e.  N=  1 for  a single-site  station.  In  general,  the 
3-N  signal  traces  in  the  edited  event  signal  were  the  vertical 
and  two  horizontal  (oriented  in  N-S  and  E-W  direction)  re- 
cordings. The  vertical  Rayleigh  wave  signal  was  obtained  from 
the  vertical  recording  while  the  radial  Rayleigh  wave  and  trans- 
verse Love  wave  signals  were  obtained  by  properly  rotating  the 
two  horizontal  recordings.  For  the  array  stations,  in  addition 
to  the  rotation  of  horizontal  recordings,  the  edited  event  signal 
was  beamformed. 

• The  edited  event  signal  with  proper  rotation  of  its  traces  (and 
beamforming  for  array  stations)  was  reformatted  to  form  the 
event  files  on  7 track  tape  to  drive  the  processing  programs  in 
the  PDP- 15/50  computer. 

• For  each  event  signal  in  the  event  files,  each  component  trace 
(Rayleigh  or  Love  wave)  was  analyzed  for  evidence  of  multi  - 
pathing  by  applying  a series  of  narrowband  filters  (as  discussed 
by  Turnbull  et  al.  , 1974;  Sun  and  Shaub,  1977).  At  each  period, 
group  velocities  were  determined  for  each  multipath  and  were 
compared  to  the  standard  group  velocity  curves  for  the  type  of 


H-5 


travel  path  encountered  (see  example  in  Sun,  1976).  Estimating 
the  correct  multipath  from  this  comparison,  the  amplitude  spec- 
tra were  obtained.  These  raw  observed  amplitude  spectra  were 
then  corrected  for  the  station  instrument  response  to  yield  the 
desired  observed  amplitude  spectra  for  the  estimation  of  the 
surface  wave  travel  path  attenuation.  These  raw  observed 
amplitude  spectra  with  appropriate  instrument  response  cor- 
rection were  further  corrected  for  geometrical  spreading  and 
travel  path  attenuation,  to  yield  the  corrected  observed  ampli- 
tude spectra.  This  corrected  amplitude  spectra  was  used  for 
the  examination  of  the  surface  wave  spectra  and  the  estimation 
of  the  seismic  source  parameters. 

In  the  above  described  surface  wave  data  reduction  procedure, 
the  two  steps  and  the  spectral  correction  in  the  third  step  are  carried  out  by 
batch  processing  on  the  IBM  360/44  computer.  The  third  step  involving  the 
narrowband  filtering  was  accomplished  using  the  Extended  Interactive  Seismic 
Processing  System  (ISPSE)  developed  by  Texas  Instruments  using  PDP-15/50 
interactive  graphics. 

C.  DATA  AVAILABILITY  AND  QUALITY 

According  to  the  work  statement  of  this  year,  it  was  requested 
that  Texas  Instruments  use  long-period  surface  wave  data  from  the  Seismic 
Research  Observatories  (SRO),  Abbreviated  Seismic  Research  Observatories 
(ASRO),  the  High  Gain  Long-Period  (HGLP)  stations  (the  same  as  the  previous 
VLPE),  and  existing  VELANET  array  data.  The  intention  of  including  as  many 
existing  long-period  recording  stations  as  possible  is  understandable  in  that  it 
provides  better  azimuthal  coverage  and  better  data  availability.  In  this  study, 
array  data  was  not  used  because  of  the  absence  of  useable  signals  caused  by 
the  distance  from  the  source  areas  of  interest,  and  in  some  cases  recording 
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malfunctions.  In  part,  the  loss  of  data  from  the  VELANET  arrays  was  com- 
pensated by  including  ASRO  data.  For  example,  some  of  the  KAA  and  MAJ 
ASRO  data  was  substituted  for  the  ILPA  and  KSRS  data  respectively.  This 
substitution  increased  the  azimuthal  coverage  to  a level  considered  satisfactory 
in  most  cases.  The  criteria  for  azimuthal  coverage  will  be  covered  later. 

After  checking  the  availability  of  the  surface  wave  recordings 
and  examining  the  results  of  the  processed  surface  wave  signals  from  all  avail- 
able recordings,  the  final  list  of  the  selected  events  was  obtained  as  shown 
in  Table  II- 1.  Table  II- 1 gives  the  event  name,  date,  origin  time,  location, 
and  PDE  m^,  together  with  the  availability  of  good  quality  surface  wave  re- 
cordings indicated  by  the  letter  'B',  'R',  or  'L'  under  the  'STATION  I.  D.  ' 
heading.  Stations  are  identified  by  the  station  number  (1.  2.  3. . . .).  'B'  means 

that  both  Rayleigh  and  Love  wave  recordings  are  useable,  'R'  means  that  only 
the  Rayleigh  wave  recording  is  useable,  and  'L'  means  that  only  the  Love  wave 
recording  is  useable.  The  descriptions  of  the  available  observation  stations 
are  given  in  Table  II-2,  including  station  name,  location,  and  identification 
number  or  letter.  It  can  be  seen  that  the  final  list  of  this  year's  selected  e- 
vents  consists  of  eight  events  (three  Eurasian  events,  four  EKZ  events,  and 
one  NTS  event)  instead  of  the  thirty  nine  events  in  the  preliminary  list.  The 
loss  of  events  from  the  preliminary  list  was  mainly  caused  due  to  the  following: 

• Bad  field  tapes  which  could  not  be  edited  for  the  desired  events, 

• Clipped  data, 

• Non- detectable  data, 

• Mixed  events  which  contain  more  than  one  event  in  the  desired 

signal  gate, 

• Malfunctions  which  could  indicate  any  non- seismic  occurrence 

which,  occurring  in  the  signal  gate,  would  mask  the  seismic 
data  under  analysis.  The  types  of  malfunctions  commonly  ob- 
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TABLE  II- 2 

DESCRIPTION  OF  OBSERVATION  STATIONS 


Station  I.  D. 

Station  Location 

# 

Name 

City  and  Country 

Lat. 

on 

Lon. 

°E 

1 

CTA 

Charters  Towers,  Australia 

-20.  0 

146.  2 

2 

CHG 

Chiang  Mai,  Thailand 

18.7 

98.  9 

w 

3 

FBK 

Fairbanks,  Alaska 

64.  9 

-148.  0 

Ou 

> 

4 

TLO 

Toledo,  Spain 

39.0 

-4.0 

A 

6 

KON 

Kongsberg,  Norway 

59.  7 

9.  6 

>4 

0 

hr1 

7 

OGD 

Ogdensburg,  New  Jersey 

41.  1 

-74.  6 

8 

KIP 

Kipapa,  Hawaii 

21.4 

-158.0 

9 

ALQ 

Albuquerque,  New  Mexico 

34.  9 

-106.  5 

10 

ZLP 

La  Paz,  Bolivia 

-16.  5 

-68. 1 

30 

ANMO 

Albuquerque,  New  Mexico 

34.  9 

-106.  5 

33 

CHTO 

Chiang  Mai,  Thailand 

18.7 

98.  9 

o 

35 

GUMO 

Guam,  Marianas  Island 

13.  6 

144.  9 

w 

36 

MAIO 

Mashhad,  Iran 

36.  3 

59.  5 

38 

NWAO 

Narrogin,  Western  Australia 

-32.9 

117. 24 

41 

TATO 

Taipei,  Taiwan 

25.0 

121.  5 

8 

52 

KAAO 

Kabul,  Afghanistan 

34.  5 

69.04 

C/5 

< 

53 

MAJO 

Matsushiro,  Japan 

36.  5 

138.  2 

B 

ALP  A 

Alaska 

65.  2 

■SB 

id 

h 

B 

LASA 

Montana 

46.7 

ESI 

< 

N 

NORSAR 

Norway 

' 

60.  8 

10.9 

served  on  the  SRO,  ASRO,  and  HGLP  data  were  glitches, 

uncorrectable  spikes,  data  dropouts,  and  calibration  signals. 

According  to  the  above  five  criteria,  only  about  twenty  per  cent  of  the  event- 
stations  corresponding  to  the  twenty  five  events  remained  to  be  used.  Fifteen 
per  cent  of  total  event- stations  were  lost  due  to  bad  field  tapes,  ten  per  cent 
due  to  malfunctions,  twenty  per  cent  due  to  mixed  events,  and  thirty  five  per 
cent  due  to  non-detectable  data.  In  addition,  nineteen  of  the  twenty- seven 
events  which  had  less  than  three  good  edited  surface  wave  traces  at  two  ob- 
servation stations  were  discarded.  This  last  condition  for  rejecting  an  event 
was  based  on  the  previous  result  that  the  source  parameters  of  a given  seismic 
event  could  be  satisfactorily  estimated  with  good  quality  surface  wave  record- 
ings at  two  observation  stations  in  two  consecutive  quadrants  (Turnbull,  et  al., 
1973).  From  the  data  availability  indicated  in  Table  U-l,  it  is  obvious  that 
the  SRO  and  ASRO  stations  present  more  and  better  quality  data  than  HGLP 
stations  do. 

In  order  to  get  better  results,  seven  processed  events  (two  from 
Eurasian  events,  two  from  NTS  events,  and  three  from  EKZ  events),  as  given 
in  Table  II- 1,  have  been  selected  from  previous  reports  to  enlarge  the  data 
base  of  this  task.  Table  II-3  presents  the  station  azimuthal  angles  with  respect 
to  those  event  locations  and  the  epicenter  distances  from  those  event  locations 
to  each  observation  station.  In  general,  the  station  coverage  for  the  selected 
events  is  thought  to  be  satisfactory  in  terms  of  station  azimuthal  distribution 
and  epicenter  distances  for  the  purposes  of  this  study. 


STATION  AZIMUTHAL  ANGLES  AND  EPICENTER  DISTANCES  FROM  THE  SELECTED  EVENTS 


SECTION  III 


THE  PROCEDURE  OF  SOURCE  PARAMETER  ESTIMATION  USING 
FAR-FIELD  SURFACE  WAVE  DATA 


Figure  III- 1 gives  a schematical  outline  of  the  procedure  to  ac- 
quire the  seismic  source  parameter  estimation  using  long-period  surface  wave 
data.  In  this  diagram,  the  following  procedures  are  used  after  selecting  events 
and  before  applying  the  amplitude  spectra  fitting  process.  The  observed  data 
are: 

• The  surface  wave  time  traces  were  obtained  by  editing  field  tapes 
with  satisfactory  event  signals. 

• The  edited  time  traces  were  analyzed  using  the  ISPSE  to  obtain 
the  demultipathed  surface  wave  spectra  (Sun  and  Shaub,  1977). 

• The  observed  surface  wave  amplitude  spectra  are  corrected  for 
instrument  response,  geometric  spreading,  and  travel  path  atten- 
uation. 

The  theoretical  data  are: 

• Earth  models  for  calculation  of  the  theoretical  amplitude  spectra. 

The  first  two  observed  data  listed  above,  were  previously  dis- 
cussed in  Section  II.  The  third  observed  data  and  theoretical  data  listed  above 
will  be  discussed  in  Subsections  IU-1  and  III-3.  Finally,  the  corrected  ob- 
served spectra  and  the  calculated  theoretical  spectra  are  input  to  an  amplitude 
spectral  fitting  process  to  obtain  estimates  of  the  source  parameters.  This 
will  be  discussed  in  Subsection  III- 2. 
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FIGURE  m-1 

SCHEMATIC  DIAGRAM  SHOWING  A PROCEDURE  TO  OBTAIN  THE  SEISMIC 
SOURCE  PARAMETER  ESTIMATION  USING  FAR-FIELD  SURFACE  WAVE  DATA 
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A. 


THE  CORRECTIONS  FOR  THE  INSTRUMENT  RESPONSE,  GEOMETRIC 
SPREADING,  AND  TRAVEL  PATH  ATTENUATION 


For  each  event,  the  observed  spectra  are  represented  by  spec- 
tral values  at  nine  discrete  periods,  namely  10  to  50  seconds  at  5 second  in- 
crements. These  spectra  have  been  equalized  at  all  observation  stations  in 
order  to  make  direct  comparison  possible.  Spectral  equalization  means  that 
raw  observed  surface  wave  spectra  have  been  corrected  for  the  station's  in- 
strument response,  geometric  spreading,  and  travel  path  attenuation.  Because 
these  three  corrections  are  well  documented  by  Sun  (1977),  no  detailed  discus- 
sions  of  the  theory  will  be  given  here. 

The  difficulty  in  correcting  for  the  instrument  response  is  not 
in  the  application,  but  in  obtaining  accurate  instrument  response  curves  for 
the  different  instruments  used  to  measure  the  surface  waves.  An  extreme  ex- 
ample of  the  difficulty  in  obtaining  accurate  corrections  is  the  HGLP  network. 
During  the  initial  assembly  and  installation  of  the  HGLP  (VLPE)  systems,  no 
attempt  was  made  to  standardize  response  characteristics.  As  a result,  each 
component  of  each  HGLP  station  had  a unique  response  curve  and  unique  digital 
sensitivity.  Moreover,  their  response  changes  frequently,  so  that  obtaining 
the  correct  instrument  response  for  a specific  time  is  a very  difficult  and  te- 
dious task.  Fortunately,  a standardized  response  curve  for  the  HLGP  system 
was  initiated  during  1976,  and  is  now  operating  for  the  HGLP  stations.  Their 
instrumentation  is  calibrated  daily  with  a step  function  input.  A sine-wave 
calibration  at  a period  of  30  seconds  is  applied  once  a month.  Calibration 
constants  and  input  levels  can  be  obtained  from  the  Albuquerque  Seismological 
Laboratory.  In  ASRO  and  SRO  systems,  the  response  characteristics  match 
the  standard  values  to  within  plus  or  minus  10%.  Digital  recording  sensitivity 
for  the  long-period  data  is  set  the  same  at  all  stations;  10,000  digital  counts 
per  micron  of  ground  displacement  at  a period  of  25  seconds  for  ASRO  systems 
and  5,000  for  SRO  systems.  Long-period  channels  are  calibrated  daily  at 


each  station  by  a 25  second  period  sine-wave  signal.  Calibration  constants  and 
input  levels  also  can  be  obtained  from  Albuquerque  Seismological  Laboratory. 


The  equation  of  correcting  for  geometrical  spreading  can  be 
written  as  follows: 


A 

op 


1/2 


(III-l) 


where: 


Acr  - amplitude  corrected  to  reference  point  r^ 

A^  - observed  amplitude  at  observation  station  p 

Ap  - epicentral  distance  between  source  and  p in  degree 

Re  - radius  of  earth  in  kilometers 

rQ  - reference  epicentral  distance  in  kilometers. 

According  to  this  equation,  it  is  obvious  that  all  the  observed  data  from  dif- 
ferent stations  of  the  same  event  will  be  corrected  to  a standard  distance  one 
kilometer  from  the  epicenter  along  different  travel  paths  to  reduce  the  effect 
of  geometrical  spreading. 

The  detailed  geographic  locations  of  the  Eurasian  events  and  the 
EKZ  events  are  displayed  in  Figures  HI-2  and  HI-3,  respectively.  F-'gure  IH-4 
shows  the  locations  of  the  available  observation  stations  in  relation  to  the  EKZ 
area  together  with  the  event-station  travel  paths.  Observed  Rayleigh  wave 
amplitude  spectra  were  available  for  nine  sampled  periods  (T  = 10  to  50  second, 
4T  = 5 sec)  of  the  seven  events  at  most  of  the  observation  stations;  while  for 
the  Love  wave,  the  observed  amplitude  spectra  were  incomplete  at  almost  half 
of  the  total  stations.  Figures  HI- 5 and  III- 6 exhibit  the  locations  of  the  avail- 
able observation  stations  in  relation  to  the  Eurasian  events  and  NTS  area,  to- 
gether with  the  event-station  travel  paths.  Except  for  the  travel  paths  from 


< t 

' t> 

I t 


HI -4 


owe  owe  owe  nx  owe  owe  ow  arse  awe  owe  me  one  iooc  lose  noe  n 

MILLER  MODIFIEO  MERCATOR  PROJECTION 
MAP  SCALE • 0.08000  INCHES/0E6REE  LONGITUDE 


FIGURE  HI-2 

LOCATIONS  OF  THE  SELECTED  EURASIAN  EVENTS 
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FIGURE  HI-3 

LOCATIONS  OF  THE  SELECTED  EKZ  EVENTS 


LOCATIONS  OF  THE  OBSERVATION  STATIONS  AND  TRAVEL  PATHS  ENCOUNTERED  TO  THE  EKZ 
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LOCATIONS  OF  THE  OBSERVATION  STATIONS  AND  TRAVEL  PATHS  ENCOUNTERED 

TO  THE  EURASIAN  EVENTS 


NTS  to  the  observation  stations  in  the  Pacific  Ocean  (such  as  KIP  and  CTA, 
which  are  entirely  oceanic),  the  travel  paths  are  either  entirely  continental 
or  mixed.  For  the  station,  CHG,  the  spectral  amplitudes  are  unusually  small, 
indicating  attenuation  possibly  caused  by  complicated  geological  features  in 
central  Eurasia;  (abnormally  thick  crust  approximately  60km  thickness  of  the 
Tibet  platform  and  the  Himalayan  fold  system).  Because  of  no  reliable  know- 
ledge of  energy  attenuation  along  specific  travel  paths,  Tryggvason’s  energy 
attenuation  curve  (Tryggvason,  1965)  has  been  used.  The  general  shape  of 
Tryggvason's  curve,  which  can  affect  the  estimation  of  the  focal  depth  and  the 
source  mechanism,  is  believed  to  be  reasonable.  Nevertheless,  by  using 
Tryggvas°n's  attenuation  correction  in  the  past,  the  seismic  moment  estimates 
of  numerous  Eurasian  seismic  events  have  fitted  very  well  with  the  seismic 
moment  versus  bodywave  magnitude  curve  determined  from  the  revised  W - 
model  given  by  Tsai  (1972b). 

B.  AMPLITUDE  SPECTRAL  FITTING  PROCESS 

In  this  Subsection,  we  abstracted  the  whole  amplitude  spectral 
fitting  process.  This  process  utilizes  the  modified  Tsai’s  method  for  the  es- 
timation of  seismic  source  parameters  (Tsai  and  Aki,  1970;  Tsai  and  Shen, 
1972;  Turnbull,  et  al.  1973).  The  modified  Tsai's  method  is  an  exhaustive 
searching  scheme  to  find  the  best  fit  between  the  observed  and  the  theoretically 
calculated  surface  wave  amplitude  spectra.  Originally,  the  modified  Tsai's 
method  was  developed  to  obtain  a source  parameter  estimate  of  an  earthquake 
modeled  by  a double-couple  source.  The  model  yields  five  source  parameters; 
seismic  moment  (M),  focal  depth  (h),  dip  (5),  slip  (X),  and  strike  (0)  angle. 

It  is  sufficient  to  say  that  the  spectral  values  of  Love  and  Rayleigh  waves  of  a 
given  frequency  depend  on  a total  of  five  parameters.  Two  of  the  five  para- 
meters, namely  the  dip  angle  6 and  slip  angle  X,  specify  the  orientation  of 
the  double -couple.  The  seismic  moment  M,  represents  the  strength  of  the 


m-io 


double-couple.  The  parameter  h,  indicates  the  depth  at  which  the  equivalent 
point  source  of  the  earthquake  dislocation  is  located  in  the  layered  half-space. 
For  an  explosion  event  with  associated  tectonic  strain  release,  the  theoretical 
source  model  which  is  used  here  combines  a point  explosive  source  and  a point 
double-couple  source.  Theoretical  surface  wave  spectra  are  generated  by  this 
combined  source  in  a layered  half-space  earth.  (Harkrider,  1964;  Ban-Menahem 
and  Harkrider,  1964;  Toksoz,  et  al.  , 1965).  The  point  explosive  source  models 
an  ideal  explosion  and  radiates  only  isotropic  Rayleigh  waves;  while  a point 
double-couple  source  models  the  explosion-associated  tectonic  strain  release 
and  radiates  anisotropic  Rayleigh  and  Love  waves.  Therefore,  the  resultant 
Rayleigh  wave  radiated  from  the  combined  source  is  the  vector  sum  of  the  Ray- 
leigh wave  from  the  point  explosion  source  and  that  from  the  point  double-couple 
source.  The  resultant  Love  wave  from  the  combined  source  comes  entirely 
from  the  point  double-couple  source.  The  resultant  amplitude  spectra  can  be 
expressed  as  follows: 

Y*  = |M  LR  (h,  5.X.0)  + M LR  (h)|.. 

1J  v JL  JL  X X I lj 

= M_  I LR  + F LR  I.. 

El  E x lij 


Yij  = ME  I LQE(h,8.A.*)|. 


(Ill- 2) 


where 


R L 

Y..  and  Y..  are  the  theoretical  Rayleigh  and  Love  wave  ampli- 
ij  ij 

tude  spectra  of  the  combined  source,  respectively, 
for  the  i-th  period  at  the  j-th  station, 

LR  and  LQ  are  the  theoretical  Rayleigh  and  Love  wave  spec- 
£ £ 

tra  of  the  double-couple  source,  respectively, 
for  a seismic  moment  of  one  dyne -cm, 


III- 11 


Mj,  and  are  the  moments  of  the  double-couple  source  and 

the  explosive  source,  respectively,  and 

F is  the  relative  strength  of  the  explosive  source 

with  respect  to  the  double -couple  source,  i.  e.  , 

F = M /M  . 
x E 

It  is  assumed  here  that  the  explosion-associated  tectonic  strain  release  takes 
place  in  the  vicinity  of  the  explosion,  i.  e.  , the  depths  of  the  explosion  source 
and  the  double-couple  source  are  the  same.  Therefore,  for  an  explosion  with 
tectonic  strain  release,  there  are  six  parameters  involved  in  the  theoretical 
source  model.  There  are  focal  depth  (h),  dip  angle  (8),  slip  angle  (\),  strike 
angle  (0),  seismic  moment  of  the  double-couple  source  (M  ),  and  the  relative 
strength  of  the  explosion  (F).  The  geometry  of  the  double-couple  source  which 
models  the  explosion-associated  tectonic  strain  release  is  shown  in  Figure  III- 7. 


Let  8^  (8,  A , <f>,  h,  F)  be  a member  in  a predefined  set  of  source 
parameters  S (8,  X,«A,h,  F).  The  exhaustive  search  in  the  amplitude  spectral 
fitting  process  tests  each  member  in  S and  obtains  the  member  which  yields 
the  best  fit  between  the  observed  and  the  theoretically  calculated  surface  wave 
amplitude  spectra.  The  entire  set  of  the  source  parameters  S can  be  con- 
veniently defined  by  the  upper  and  lower  bounds  of  each  source  parameter 
with  discrete  increment,  that  is: 


h,  < h < h 
1 m 

8.  - 6 - 8. 

1 k 

Xj  < X S X^ 

<f>  < ♦ < 6 

1 n 

F.  < F < F 

i q 


with  increment  of  Ah 
with  increment  of  A8 
with  increment  of  AX 
with  increment  of  A0 
with  increment  of  AF 


where  m,  k,  / , n,  and  q are  positive  integers.  Thus,  the  entire  set  S is  a 
discrete  set  consisting  of  m«k*/*ri*q  members.  In  order  to  cut  down  com- 


8 = Dip  angle 

A = Slip  angle 

^ = Strike  angle 

0 = Azimuth  angle  with  respect  to  strike  direction 

n = Surface  normal 
A,^=  Orthogonal  vectors  on  the 
surface  of  fault  plane 
R = One  force  vector  of  double- 
couple force  on  the  surface 
of  fault  plane 

0 = Arbitrary  observation  station 

TT  = Geographical  north 

FIGURE  III -7 

GEOMETRY  OF  THE  DOUBLE-COUPLE  SOURCE  WHICH  MODELS 
THE  EXPLOSION -ASSOCIATED  TECTONIC  STRAIN  RELEASE 
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puter  expenses  and  widen  the  searching  range  of  some  parameters,  a binary 
method  was  applied  to  modify  the  amplitude  spectral  fitting  process.  Accord- 
ing to  previous  experience,  the  reasonable  range  and  increment  of  each  para- 
meter were  used  as  follows: 

30°  < 5^90°,  Ah  - 10°;  - 90°  < A < 90°,  Ak  ~ 30°; 

0°  <4  -175°,  At  = 5°;  0 $ F < 3.  0,  AF  = 0.  3; 
h depends  on  the  different  earth  models  (see  Table  III-l). 

The  non-uniform  increment  in  h is  dictated  by  the  individual 
layer  thickness  of  the  layered  half-space  earth  model  which  is  chosen  for  the 
source  region.  Three  parameters,  which  are  depth,  strike  angle  and  F value, 
were  involved  in  the  binary  method.  Therefore,  we  can  enlarge  the  above  in- 
crements by  2n  times  in  which  n indicates  the  desired  number  of  passes. 

The  same  search  strategy  was  applied  to  all  of  the  events  analyzed;  both  earth- 
quakes and  explosions,  were  searched  through  the  same  kind  of  range  and  in- 
crement of  each  parameter  to  get  final  minimum  error  parameter  estimations. 

R L 

Let  X..  and  X..  denote  the  observed  Rayleigh  and  Love  wave 
ij  iJ 

spectral  amplitudes,  respectively,  for  the  i-th  period  at  the  j-th  station. 

The  fitting  error  taken  as  the  difference  between  the  observed  spectra  and  the 
corresponding  theoretical  spectra  is  referred  to  as  the  residual  e,  defined  as: 
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Substituting  equation  (III- 2)  into  equation  (ni-3),  we  have: 
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TABLE  m-1 

RANGES  AND  INCREMENTS  OF  THE  SEARCHED  DEPTH 


N*Earth 

>Model 

DepthV 


EKZ 

B-C  Foldbelt 


4.0 

6.0 

8.0 

9.  5 

11.5 

14.0 

16.0 

18.0 

20.0 

21.  5 

25.0 

30.0 

35.0 

40.0 

43.0 

45.0 
47.  5 
49.  5 


Eurasian 

G-B 

NTS 

35-CM 

0.  5 

0.  5 

2.0 

1.  5 

4.0 

2.  5 

6.0 

4.0 

8.0 

6.  0 

10.0 

8.0 

12.0 

10.0 

14.0 

12.0 

16.0 

14.0 

18.0 

16.0 

19.5 

18.0 

20.  25 

20.0 

22.  5 

22.  5 

25.0 

25.0 

27.  25 

27.  5 

30.0 

30.0 

32.  5 

32.  5 

35.0 

34.  5 

37.  5 

37.  5 

39.  5 

42.  5 

42.0 

47.0 

45.0 

50.0 

47.  5 

49.  5 

* Earth  Model:  details  discussed  in  Section  IH-C 


In  the  amplitude  spectra  fitting  process,  first  the  moment  is  estimated  by 
lizing  the  residual  «(s^forall  i and  j,  or  explicitly: 


minimi  j 


M. 


z z 

[x* 

LRe  + F LRx  | 

.,+x1; 
ij  ij 

lqe 

■ 

ij 

z z 

j i 

’LRI 

+ F LR 

X 

2.  + 
ij 

|lqe 

2 

ij 

3 

(in- 5) 


In  other  words,  for  a particular  member  s^  in  the  set  of  source  parameters 
S,  the  optimum  moment  is  first  found  from  equation  (III- 5).  Then,  the  result- 
ant M_  (s  ) is  used  in  equation  (IU-4)  to  compute  the  minimum  residual  (s  ) 
h,  p p 

for  this  particular  member  s^.  The  minimum  residual  €(8^)  is  computed  for 
each  of  the  members  in  S,  This  e value  is  used  as  a convenient  criterion  for 
judging  the  goodness  of  the  fit  between  the  observed  and  the  theoretically  cal- 
culated spectra  for  each  member  s^.  In  the  end,  the  source  parameter  es- 
timation of  a given  event  will  be  the  member  s0^”14^  whose  theoretical 

P 

surface  wave  spectra  fits  moat  closely  the  observed  spectra ; that  is: 


€ min  ^ s°ptimal)  = Minimum  [e(»p)]  for  all  • in  S.  (HI-6) 

In  order  to  attain  a greater  degree  of  confidence  in  using  the 
source  parameter  estimation  based  on  the  minimum  residual  criterion,  a de- 
tailed knowledge  of  the  residual  distribution  about  the  absolute  minimum  given 
by  equation  (III-6)  would  be  helpful  to  guard  against  a possible  spurious  mini- 
mum generated  by  poor  quality  data.  Therefore,  in  the  amplitude  spectral 
fitting  process  used  here,  distributions  of  <(*p)  with  respect  to  the  source 
parameters  are  examined. 

In  the  fitting  process,  the  focal  depth  is  varied  last,  because 
the  focal  depth  (h)  will  have  the  most  affect  on  the  shape  of  the  theoretical  sur- 
face wave  amplitude  spectra,  especially  for  Rayleigh  waves  (Turnbull,  et  al., 
1974).  There  is  one  €min  at  one  focal  depth  h^,  that  is: 
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h. 

6 1 
min 


= Minimum  <<[s  (6,k,<P,  h.,  F)ll  for  all  6,  \,4>  and  F 
in  ( P l f 


at  depth  h^.  (HI -7) 

The  variation  of  with  respect  to  the  focal  depth  is  obtained  and  a standard 

fh  * hi  • 

! \ I ^ . ,,  is  determined  from  those  f \ for  the  purpose 

minj  standard  min  K v 

of  picking  appropriate  members  in  the  set  of  source  parameters  S to  form  a 

sample  set  Sd.  In  the  entire  set  of  source  parameters  S,  those  members 

whose  e(s  )'s  are  less  than  [e  * 1 , , are  included  in  the  sample  set 

p L rmi\J  standard  r 

Sd.  These  form  the  distributions  of  c(s^)  with  respect  to  the  source  para- 
meters. These  residual  distributions  for  the  various  source  parameters  will 
be  more  easily  understood  from  the  results  which  will  be  presented  later  for 
the  source  parameter  estimations  of  the  selected  events. 


C.  EARTH  MODEL 

From  the  theoretical  part  of  Figure  III- 1 , it  is  obvious  that  a 
layered  half-space  earth  model  (which  describes  the  source  region  earth  struc- 
ture) has  to  be  obtained  first  in  order  to  generate  theoretical  spectra.  Then, 
we  can  use  earth  models  to  calculate  theoretical  surface  wave  amplitude  spec- 
tra which  are  needed  in  the  amplitude  spectral  fitting  process.  In  this  task, 
three  earth  models,  as  shown  in  Figure  III-8,  are  selected  for  the  EKZ,  NTS, 
and  Eurasian  events. 

Both  groups  of  EKZ  events  and  NTS  events  employed  the  same 
earth  models  that  were  used  in  the  previous  report  (Sun,  1977).  For  the  EKZ 
events,  the  earth  model  is  constructed  from  a representative  subsurface 
crustal  structure  of  the  Balkhash -Chingiz  Foldbelt,  (see  Figure  III-8a).  Since 
these  crustal  cross  sections  only  extend  to  a depth  of  approximately  50  km, 
this  model  was  created  by  merging  the  crustal  cross  section  with  the  Gutenberg- 
Bullen  model  (Turnbull,  et  al.,  1974)  at  a depth  of  50  km.  For  the  NTS  events, 
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the  combined  earth  model  (which  is  the  model  that  merges  the  35CM2  models 
of  Alexander  (1963)  with  the  Gutenberg-Bullen  model  (see  Figure  III-8c)  at  a 
depth  of  150  km)  was  utilized.  Actually,  the  35CM2  model  (see  Figure  IH-8b) 
is  a simple  structure  consisting  of  about  2-3  km  of  'sedimentary'  material 
overlying  a rather  uniform  section  of  'granite'  material  approximately  23  km 
thick.  The  section  of  granite  in  turn  overlies  an  intermediate  section  approx- 
imately 25  km  thick.  For  the  Eurasian  events,  it  was  assumed  that  the  normal 
Gutenberg-Bullen  earth  model  is  most  appropriate  (Turnbull,  1974). 
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SECTION  IV 


RESULTS  AND  CONCLUSIONS  OF  SOURCE  PARAMETER 
ESTIMATIONS  FOR  EKZ,  EURASIAN,  AND  NTS  EVENTS 


In  the  previous  section,  we  discussed  the  procedures  to  obtain 
the  theoretical  amplitude  spectra  and  the  observed  amplitude  spectra  from 
earth  models  and  measured  seismic  waveforms.  Also,  we  described  the  mod- 
ified amplitude  spectra  fitting  process  which  is  utilized  to  get  seismic  source 
estimations.  In  this  section,  the  results  of  the  seismic  source  parameter  esti- 
mation will  be  presented  separately  for  EKZ,  Eurasian  and  NTS  events.  The 
results  of  using  these  source  parameter  estimates  for  discrimination  between 
earthquakes  and  explosions  will  be  examined.  Finally,  a brief  discussion  of 
all  the  results  also  will  be  given. 

A.  THE  EKZ  EVENTS 

For  EKZ  events,  observed  surface  wave  data  are  available  for 
both  Rayleigh  and  Love  waves  at  most  stations,  except  for  the  event  EKZ/427/ 

75  in  which  only  Rayleigh  wave  data  are  used  at  all  available  stations.  However, 
Sun  (1976)  compared  the  results  obtained  from  spectral  fittings  using  both  Ray- 
leigh and  Love  waves  to  spectral  fittings  using  Rayleigh  waves  only.  He  re- 
ported no  significant  difference  in  the  source  parameter  estimations.  There- 
fore, the  event  EKZ/427/ 75  is  still  in  our  final  event  list. 

The  optimal  solution  of  source  parameters  based  on  the  minimum- 
residual  criterion  are  listed  in  Table  IV- 1.  The  residual  distributions  with  re- 
spect to  the  estimated  depth,  dip  angle,  slip  angle,  and  strike  angle  are  given 
in  Figure  IV-1  through  Figure  IV-7  with  the  spectral  fits  for  each  event.  As 
described  in  more  detail  by  Sun  (1976),  the  residual  distribution  for  one  particular 
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source  parameter,  for  example,  depth,  consists  of  the  variation  of  the  mini- 
mum-residual values  with  respect  to  the  depths  tested  in  the  spectra  fitting 
process  and  the  percentage  of  occurrence  of  minimum-residual  values  at  the 
various  tested  depths.  For  each  event,  the  residual  distributions  are  given 
in  part  (a)  of  Figure  IV- 1 through  Figure  IV-7.  The  best  spectral  fits  of  Ray- 
leigh and  Love  waves  (as  defined  by  equation  (III- 6)^  are  shown  in  part  (b)  of 
Figure  IV-1  through  Figure  IV-7.  In  these  figures,  the  residual  values  are 
normalized  with  respect  to  the  residual  error  of  the  optimal  minimum  var- 
iance solution  which  most  closely  fits  the  observed  spectra.  These  are  indi- 
cated as  the  relative  minimum  residuals  (R.  M.  R.) . The  observed  spectra 
are  indicated  by  crosses  (x)  and  the  best  fit  theoretical  spectra  are  represent- 
ed by  open  circles  (o)  connected  by  straight  lines.  When  there  are  no  observed 
spectra  for  either  the  Rayleigh  or  Love  wave,  the  crosses  are  plotted  on  top 
of  the  corresponding  open  circles.  From  the  residual  distributions  shown  in 
part  (a)  of  Figures  IV-1  through  Figure  IV-7,  the  estimates  of  the  source 
parameters  of  each  event  can  be  represented  in  terms  of  the  probable  range  of 
the  estimate  and  the  percentage  frequency  of  that  probable  range.  The  estima- 
tions of  source  parameters  based  on  the  residual  distributions  are  listed  in 
Table  IV- 2.  To  give  a clearer  explanation  of  Table  IV-2,  consider  the  event 
EKZ/03/20/76.  It  has  a depth  estimate  of  8.  0 km  indicated  in  Table  IV-1; 
while  Table  IV-2  indicates  that  100%  of  the  members  in  the  sample  set  s^  of  the 
tested  source  parameters  (as  described  in  Subsection  III- 3)  indicate  that  the 
probable  range  of  the  depth  estimate  is  from  8.  0 km  to  11.5  km. 

Referring  to  Table  IV-1,  these  seven  events  can  be  divided  into 
two  groups  according  to  the  dip  and  slip  angle  estimates  of  the  tectonic  strain 
release  associated  with  these  events.  The  first  group  is  nearly  dip- slip  for 
one  event  and  purely  dip- slip  for  four  events,  and  the  second  group  is  nearly 
strike-slip  for  two  events.  These  results  correspond  to  those  of  Sun  (1977). 

The  estimated  dip  angles  are  quite  consistent  from  event  to  event,  except  for 
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the  event  EKZ/ 723/73.  The  reason  for  this  variation  may  be  caused  by  the 

poor  spectral  fits  for  the  Rayleigh  wave  data  of  LASA  station.  The  probable 

o o 

range  of  dip  angles,  as  given  in  Table  IV- 2,  are  approximately  40  to  60  . 

For  the  estimated  strike  angles,  there  is  no  consistency  among  the  seven  e- 

vents.  The  majority  of  the  depth  estimates  of  these  events  are  shallow;  from 

0.  5 km  to  4.  0 km  except  8.  0 km  for  both  events  EKZ/03/20/76  and  EKZ/12 

14/  3 . The  shallow  depths  are  coincident  with  the  fact  that  those  events  are 

presumed  underground  explosions.  For  the  event  EKZ/03/20/76,  (one  of  the 

two  events  of  8.  0 km  depth),  Strauss  (1977)  reported  that  the  values  for 

this  event  as  recorded  at  ANMO,  GUMO,  MAZO,  and  NWAO  all  fall  in  the 

earthquake  population,  indicating  this  event  is  an  earthquake.  In  our  M ver- 

£ 

sus  m plot,  which  is  discussed  later,  this  event  also  falls  in  the  earthquake 
population.  These  descriptions  might  explain  the  8.  0 km  depth  of  this  event. 
The  probable  reason  for  the  other  anomalous  event  EKZ/ 1214/3  will  be  dis- 
cussed in  Subsection  IV-E.  The  F-values  vary  from  0.  0 to  0.  6 among  these 
seven  events.  This  indicates  that  the  double-couple  source  is  for  the  most 
part  responsible  for  the  observed  surface  waves  of  these  events.  On  physical 
grounds,  this  might  be  expected  due  to  the  expected  order  of  magnitude  drop 
in  the  amplitude  of  compressional  waves  emitted  from  a shallow  explosion 
source  and  possibly  from  such  a shallow  strain  release  source  due  to  the  free 
surface  reflection. 

Judging  the  spectral  fits  visually,  these  seven  events  are  fairly 
good  at  most  available  stations.  In  general,  the  Rayleigh  wave  amplitude  spec- 
tral fits  are  better  than  those  of  the  Love  wave. 

B.  THE  EURASIAN  EVENTS 

For  the  selected  Eurasian  events,  the  estimations  of  source 
parameters  based  on  the  minimum- residual  criterion  are  listed  in  Table  IV-3. 
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Those  based  on  the  residual  minima  frequency  of  occurrence  distributions  are 

listed  in  Table  IV-4.  The  residual  distributions  with  respect  to  the  estimated 

depth,  dip  angle,  slip  angle,  and  strike  angle  are  given  in  Figure  1V-8  through 

Figure  IV- 12  with  the  spectral  fit  for  each  event.  i 

Referring  to  Table  IV-3,  these  five  Eurasian  events  may  be  di- 
vided into  two  groups  according  to  their  focal  depth.  In  group  1,  ECU/02/03/ 

72  and  TZK/03/17/72  have  the  deeper  focal  depth  in  the  range  of  20  to  45  km 
and  in  group  2,  the  other  three  have  shallow  focal  depths  of  approximately  4.  0 
km  to  8.  0 km.  In  the  first  group,  it  is  interesting  to  note  the  larger  F-values 
of  the  two  deeper  events.  The  larger  F-values  may  suggest  that  more  energy  1 

is  released  from  the  dilatational  point  source  than  from  that  of  the  double- 
couple point  source  in  the  deeper  earthquakes.  In  the  second  group,  event  NVZ/ 

09/01/77  would  be  classified  as  an  explosion  due  to  its  special  origin  time  and 

its  seismic  moment  falling  into  the  explosion  population  of  our  M versus  m 

Ed 

plot. 

For  event  ECU/02/03/72,  the  spectral  fits  yielded  a focal  depth 
in  the  range  of  42  to  47.  5 km  with  the  best  fit  at  45  km.  Comparing  this  result 
with  the  depth  of  39  km  indicated  in  the  NEIS  list,  the  difference  is  not  signifi- 
cant. The  best  fit  dip  and  slip  angles  suggested  an  oblique- slip  source  mechan- 
ism with  a strike  direction  of  N15°E.  The  estimated  F-value  is  0.75. 

For  event  TZK/03/17/72,  the  optimal  solution  showed  a dip-slip 
fault  mechanism  with  an  approximate  strike  direction  of  N55°E.  The  estimated 

' 

focal  depth  lies  in  the  range  of  18  to  22.  5 km  with  the  best  fit  of  depth  at  19.  5 

r 

km.  The  difference  between  this  estimated  depth  and  the  depth  of  26  km  (indi- 

r 

cated  in  the  NEIS  list  ) might  be  caused  by  the  poor  quality  of  the  amplitude 
spectral  fit  at  10- second  periods,  because  the  variation  of  the  amplitude  spec- 
tra at  shorter  periods  (less  than  25  seconds)  will  affect  the  estimation  of  focal 
depth  mostly.  The  estimated  F-value  is  1.6. 
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For  event  LBR/08/10/77,  the  focal  depth  is  located  in  the  range 
of  2.  0 to  4.  0 km  with  the  greatest  possibility  at  4.  Okm.  The  focal  depth  obtained 
from  bodywave  depth  phases  is  2.  0 km.  Only  two  observation  stations  were 
available  for  this  event;  the  best  fit  theoretical  spectra  followed  quite  closely 
the  general  shapes  of  their  observed  counterparts.  And  this  yielded  an  oblique- 
slip  fault  mechanism  with  a strike  direction  of  N120°E.  The  F-value  is  0.  6. 

For  event  UZB/07/14/77,  8.0  km  depth  is  given  in  Table  IV-3. 
The  spectral  fits  suggested  a dip- slip  mechanism  with  a strike  angle  of  N40°E. 

For  event  NVZ /09/01 /77,  as  discussed  previously,  this  event 

might  be  an  explosion  from  Novaya  Zemlya.  The  estimated  focal  depth  is  4.  0 

km.  The  spectral  fits  yielded  a dip- slip  fault  mechanism  with  a strike  direc- 
o 

tion  of  N15  E.  The  estimated  F-value  is  0.  1. 

C.  THE  NTS  EVENTS 

For  the  selected  NTS  events,  the  estimations  of  source  param- 
eters based  on  the  minimum- residual  criterion  are  listed  in  Table  IV- 5. 

Those  based  on  the  residual  minima  frequency  of  occurrence  distribution  are 
listed  in  Table  IV- 6.  The  residual  distributions  with  respect  to  estimated 
depth,  dip  angle,  slip  angle,  and  strike  angle  are  given  in  Figure  IV- 13  to  IV- 
15  with  the  spectral  fit  for  each  event. 

According  to  Table  IV-5,  a dip-slip  fault  mechanism  was  sug- 
gested for  these  three  events.  The  estimated  focal  depths  are  very  shallow 
at  1.  5 km  for  events  NT1/317/76  and  NTS/02/07/76.  The  strike  angle  esti- 
mates are  quite  close  to  the  N-S  direction  as  assumed  by  some  authors  (Toksoz, 
et  al. , 1965;  Lambert,  et  al. , 1972;  Sun,  1977).  From  the  estimated  seismic 
moment,  the  event  NTS/02/07/76  would  be  classified  as  an  earthquake,  and  the 
reason  for  this  is  twofold.  First,  this  event  fell  into  the  earthquake  population 
of  our  Mg  versus  m^  plot,  and  secondly,  the  4.  5 km  depth  is  indicated  in  the 
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NEIS  list.  Events  NTS/212/76  and  NT1/317/76  are  the  presumed  explo- 
sions. The  estimated  F-values  varied  from  0.  5 to  1.  25. 

D.  SEISMIC  MOMENT  VERSUS  mL  PLOTS 

D 

As  we  mentioned  in  Section  III,  the  seismic  moment  is  used  as 
a parameter  for  minimizing  the  error  in  our  spectral  fitting  procedure.  Thus, 
the  seismic  moment  is  determined  along  with  other  source  parameters  in  the 
process.  Figures  IV-16  (a),  (b),  and  (c),  display  the  estimated  seismic  mom- 
ent as  a function  of  the  NEIS  mfa  magnitude  for  each  of  the  different  groups  of 

events.  In  these  figures,  the  solid  line  L_  represents  the  theoretically  pre- 
^ £ 

dieted  relation  between  the  seismic  moment  and  the  event  m.  . This  relation 

2 b 

was  based  on  the  revised  CJ  - model  determined  by  Tsai  (1972a)  for  earth- 
quake source  spectra.  Another  solid  line  L^.  for  the  explosion  population  is 
a straight  line  parallel  to  the  theoretically  predicted  line  L for  the  earth- 
quake. 

In  Figure  IV-16  (a),  we  used  M , which  was  the  seismic  mom- 
ent  obtained  from  the  spectra  fitting  process  using  the  combined  source  for 
each  event,  for  both  the  earthquake  population  and  explosion  population.  Ac- 
cording to  event  characteristics,  the  estimated  moments  of  all  events  except 
for  the  presumed  explosion  EKZ /06/29/77,  fell  in  two  well  separated  popula- 
tions and  correlated  closely  with  the  two  theoretically  predicted  curves.  In 
Figure  IV-16  (b),  the  total  moment,  M^,  which  combined  the  moment  of  the 
double-couple  source  and  the  moment  of  the  explosive  source,  was  used  in 
stead  of  M_  in  Figure  IV-16  (a).  It  is  noticeable  that  the  total  moment  gives 
better  discrimination  than  the  moment  of  double-couple  source  only.  In  Fig- 
ure IV-16  (c),  the  previously  obtained  seismic  moment  estimates  of  various 
central  Eurasian  earthquakes  (Tsai  and  Shen,  1972;  Turnbull,  et  al. , 1973, 

1974;  Sun,  1976)  and  the  selected  NTS,  EKZ  and  PNE  events  (Sun,  1977)  were 
merged  into  Figure  IV-16  (a)  to  yield  this  plot.  The  legend  on  the  left  corner 
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of  Figure  IV-16  (c)  only  described  the  merged  data.  The  present  data  applied 
in  this  plot  still  use  the  same  symbols  as  shown  in  Figure  IV-16  (a). 

Based  on  examining  these  three  plots,  the  following  remarks 

can  be  stated: 

• The  total  moment  (M^,)  seems  to  offer  a better  discrimination. 

• The  separation  between  earthquake  and  explosion  populations  is 

approximately  one  m,  unit.  This  coincides  with  the  previous 

b 

results  (Sun  1977). 

• Less  scattering  in  both  populations  keeps  these  plots  good  dis- 
criminants. 

E.  CONCLUSIONS 

The  following  remarks  summarize  the  results  in  the  preceding 
subsections  concerning  source  parameter  estimates  obtained  from  surface 
wave  data  using  spectral  fitting  procedures: 

• Compared  with  the  previous  results  (Tsai  and  Shen,  1972;  Turn- 
bull  et  al. , 1973,1974;  Sun,  1976,1977),  the  spectral  fitting  of  this 
study  offered  a very  promising  result.  Even  the  LQ  spectra  at 
most  available  stations  produced  an  acceptable  spectral  fit  for 
most  of  the  events.  For  most  events,  the  average  spectral  am- 
plitude indicated  the  same  aspect,  namely,  higher  at  shorter  per- 
iods and  decreasing  toward  longer  periods.  According  to  the  ex- 
amination of  the  theoretical  spectral  variation  with  respect  to 
focal  depth  for  the  combined  source,  this  aspect  suggested  a very 
shallow  focal  depth  (Turnbull  et  al. , 1974).  This  speculation  was 
proven  by  results  of  depth  estimates  in  Tables  IV-1,  IV-3,  and  IV-5. 

• Figure  IV-1 7 shows  the  percentage  of  events  occurring  at  each 
estimated  slip  angle  of  the  double-couple  source  of  tectonic 
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strain  release.  All  three  of  the  NTS  events  were  indicated  as 
dip  slip.  For  the  EKZ  events  and  Eurasian  events,  the  results 
are  mixed  with  no  difference  of  mechanism  discernable  between 
both  sets  of  events.  For  Eurasian  events,  obviously  no  general 
trend  among  those  slip  angles  was  expected  due  to  these  events 
being  geographically  scattered  around  the  Eurasian  continent. 

In  general,  the  estimates  of  focal  depth  of  the  selected  earth- 
quakes closely  agree  with  those  determined  from  bodywave  depth 
phases.  Figure  IV-18  shows  the  percentage  of  all  of  the  presum- 
ed explosions  occurring  at  each  depth.  This  distribution  indi- 
cates a noticeable  rightward  skew  toward  deeper  depth  estimates. 
The  mean  value  of  this  distribution  is  located  at  a 3.5  km  depth. 
Most  events  occur  in  the  confidence  area  of  this  distribution; 
that  is  an  area  of  depth  occurring  at  1.  5 km  to  4.  0 km.  For 
events  NTS/212/76  and  EKZ/1214/3,  the  estimated  depths 
are  6,  0 km  and  8.  0 km,  respectively.  The  reason  for  this  might 
be  a close  association  with  irregularities  in  crustal  structure 
and  hence,  with  anomalous  crustal  rigidity  or  strength  right  be- 
neath the  shot  point. 

Figure  IV-19  shows  the  percentage  of  total  events,  total  explo- 
sions, and  EKZ  events  occurring  at  each  dip  angle.  The  mean 
values  and  standard  deviations  of  these  three  distributions  are 
shown  in  Table  IV-7.  According  to  Table  IV-7  and  Figure  IV- 
19,  the  estimates  of  dip  angle  of  nine  explosions  most  likely 
occur  between  40°  to  50°.  This  kind  of  consistency  might  imply 
that  the  explosion  can  induce  structure  rupture  along  the  planes 
of  maximum  shear  stress  which  usually  are  oriented  around  45°. 

Figure  IV- 20  shows  the  percentage  of  explosions  and  earthquakes 
at  each  F-value.  No  consistency  between  the  explosion  population 
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THE  FREQUENCY  DISTRIBUTIONS  OF 
EXPLOSIONS  AND  EARTHQUAKES  AT  EACH  F-VALUE 


and  the  earthquake  population  has  been  found.  In  general,  the 
estimates  of  F-value  of  explosions  are  smaller  than  those  of 
earthquakes.  This  result  is  quite  different  from  what  we  ex- 
pected. The  large  F-value  of  earthquakes,  especially  for  the 
two  deeper  events  ECU/02/03/72  and  TZK/03/1 7/72,  might  be 
caused  by  the  volume  change  during  the  rupture  at  deep  focal 
depth.  This  change  usually  favors  an  explosive  point  source. 

The  small  F-value  of  explosions  means  that  the  double-couple 
source  is  dominant  at  the  periods  considered.  The  association 
of  higher  F-values  with  deeper  earthquakes  suggests  that  a high- 
er dilatational  component  is  required  to  overcome  the  overbur- 
den pressure.  The  lack  of  dilatation  of  explosions,  random 
strike  and  deeper  focus  suggests  that  the  double-couple  source 
of  explosions  may  be  associated  with  large  dislocations  on  a 
fracture  zone  surrounding  the  source.  Radiation  from  the  deep- 
er part  of  the  fracture  zone  is  more  visible  due  to  the  effect  of 
the  free  surface.  Compared  with  the  previous  results  (Sun, 
1977),  less  consistency  and  small  values  of  the  estimates  of  F- 
value  have  been  found  in  this  study.  The  reason  for  this  dif- 
ference is  that  we  used  the  binary  method  to  search  through 
wider  ranges  and  finer  increments  of  parameters  such  as  depth, 
F-value,  and  strike  angle,  without  any  constraint.  Also,  Sun's 
results  were  obtained  by  apparently  limiting  the  search  for 
minimum  residual  F-values  to  values  greater  than  one.  This 
explains  his  numerous  occurrences  of  F=1  as  edge  phenomena 
in  the  exhaustive  search. 
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SECTION  V 
RECOMMENDATIONS 


It  is  recommended  that  future  studies  be  directed  toward  the 
following  points: 

• Confirm  the  apparently  universal  surface  wave  moment  versus 
m^  relationship  with  a much  larger  data  base  of  earthquakes 
and  explosions. 

• Evaluate  surface  wave  moment  as  a scaling  parameter  for  esti- 
mating yields  of  explosions. 

• Evaluate  the  sensitivity  of  source  mechanism  solutions  to  the 
distribution  of  stations  and  path  corrections. 

• Check  the  consistency  of  combined  source  mechanisms  on  a 
large  earthquake  and  aftershock  sequence. 

• Develop  feasible  procedures  for  applying  the  combined  source 
mechanism  as  an  extended  event  identification  technique. 
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